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Abstract 

This paper reports the cross section measurements for the process ep — > e J/tp p for 
Q 2 < 4 GeV 2 at = 296 GeV, based on an integrated luminosity of about 0.5 pb _1 , 
using the ZEUS detector. The J/ip was detected in its e + e~ and decay modes. The 

photoproduction cross section was measured to be 52 ^ 10 n b at an average 7p centre 
of mass energy of 67 GeV and 71^0 i 12 nb at 114 GeV. The significant rise of the cross 
section compared to lower energy measurements is not in agreement with VDM models, 
but can be described by QCD inspired models if a rise in the gluon momentum density 
at low x in the proton is assumed. 
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1 Introduction 



Elastic photoproduction of J/ip (7P^ J/V'P) is particularly interesting as the production cross 
section can be calculated as a function of the 7p centre of mass (cm.) energy, W, both using 
the Vector Dominance Model[[L| (VDM) extended to the heavy quarks 0, |3|, |J, and, because 
of the large value of the J/ip mass (Mj/^,), with QCD inspired models || ||. Figures |T|sl and 
|TJb show the elastic J/ip photoproduction mechanisms according to VDM and QCD inspired 
models, respectively. A characteristic of these QCD inspired models is that the cross section 
is proportional to the square of the proton's gluon momentum density. The typical x range 
probed here, where x (~ M 2 ^/W 2 ) is the fraction of the proton momentum carried by the 
gluon, is approximately 5 x 10~ 3 to 5 x 10~ 4 , corresponding to a W range between 40 and 140 
GeV. For gluon distributions increasing at low x, the QCD approaches predict much higher 
cross sections than those from the VDM in this W region. 

Previous measurements of the photoproduction of J/ip are available in the W range between 4 
GeV and 28 GeV|7], ||. However, some results include elastic as well as other production mecha- 
nisms of the J/i/j. The HI experiment has also reported|| a J/ip photoproduction cross section 
measurement at an average cm. energy of 90 GeV, that contains an unknown contribution 
from inelastic J ftp photoproduction. 

This paper reports the measurements of photoproduction cross sections of elastically produced 
J/^'s in the reaction ep — > e J/ip P, followed by J/V> — > e + e~or J/i/; — > fi + fi~, with the ZEUS 
detector. The present data sample contains events with Q 2 < 4 GeV 2 , with a median Q 2 ~ 
10~ 3 GeV 2 , in the W range between 40 and 140 GeV. Neither the scattered electron nor the 
proton is detected in this measurement. The contribution of J/ip production where the proton 
diffractively dissociated is subtracted to obtain the elastic photoproduction cross section. 



2 HERA 

The data presented were collected during the 1993 running period of HERA and represent 
an integrated luminosity of about 0.5 pb _1 . HERA operated in 1993 with an electron beam 
energy of 26.7 GeV and a proton beam energy of 820 GeV. A total of 84 colliding electron 
and proton bunches was used; in addition 10 electron and 6 proton unpaired bunches were 
used for background studies. The time between bunch crossings at HERA is 96 ns. A typical 
instantaneous luminosity of ~ 6 x 10 29 cm~ 2 s _1 was delivered. 



3 The ZEUS detector 



The ZEUS detector[]T(J is a hermetic, general purpose magnetic detector with a tracking region 
surrounded by a high resolution calorimeter followed in turn by a backing calorimeter and the 
muon detector. A short description of the components relevant to this analysis is given here. 
They are the vertex detector [11], the central-tracking detector [12], the uranium-scintillator 
calorimeter ||13|| and the barrel and rear muon detectors ||14j| . 

Charged particles are measured by the ZEUS inner tracking detectors, which operate in a 
magnetic field of 1.43 T provided by a thin superconducting coil. Immediately surrounding the 
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beampipe is the vertex detector (VXD) consisting of 120 radial cells, each with 12 sense wires. 
It uses a slow drift velocity gas and the presently achieved resolution in the XYf\ plane is 50 fim 
in the central region of a cell and 150 jum near the edges. Surrounding the VXD is the central 
tracking detector (CTD) which consists of 72 cylindrical drift chamber layers, arranged in 9 
'superlayers'. With our present calibration of the chamber, the resolution of the CTD is around 
260 /im. The resolution in transverse momentum for tracks going through all superlayers is 
v{Pt)/pt ~ y (0.005) 2 p^ + (0.016) 2 where pt is in GeV. The single hit efficiency is greater 
than 95%. The efficiency for assigning hits to tracks depends on several factors: very low pt 
tracks suffer large systematic effects which reduce the probability of hits being assigned to them, 
and the 45° inclination of the drift cells also introduces an asymmetry between positive and 
negative tracks. Nevertheless, the track reconstruction efficiency for tracks with p? > 0.1 GeV 
is greater than 95%. Using the combined data from the VXD and CTD, resolutions of 0.4 cm 
in Z and 0.1 cm in radius in the XY plane are obtained for the primary vertex reconstruction. 
From Gaussian fits to the Z vertex distribution, the r.m.s. spread is found to be 10.5 cm, in 
agreement with the expectation based on the HERA proton bunch length. 

The high resolution uranium-scintillator calorimeter (CAL) covers the polar angle range be- 
tween 2.2° < 9 < 176.5°, where 9 = 0° is the proton beam direction. It consists of three parts: 
the rear calorimeter (RCAL), covering the backward pseudorapidity^ range (—3.4 < rj < —0.75); 
the barrel calorimeter (BCAL) covering the central region (—0.75 < i] < 1.1); and the forward 
calorimeter (FCAL) covering the forward region (1.1 < rj < 3.8). The calorimeter parts are 
subdivided into towers which in turn are subdivided longitudinally into electromagnetic (EMC) 
and hadronic (HAC) sections. The sections are subdivided into cells, each of which is viewed 
by two photomultiplier tubes. Under test beam conditions the CAL has an energy resolution, 
in units of GeV, of oe = 0.35y / £'(GeV) for hadrons and <je = Q.18\J E(GeV) for electrons. The 
CAL also provides a time resolution of better than 1 ns for energy deposits greater than 4.5 
GeV, and this timing is used for background rejection. 

The muon detectors, placed outside the calorimeter, are also divided into three sections, covering 
the forward, barrel and rear regions. In the barrel and rear regions limited streamer tube (LST) 
chambers before (inner) and after (outer) an 80 cm thick magnetized iron yoke are used. Only 
the inner chambers of the barrel and rear muon detectors (BMUI and RMUI) were used for the 
present analysis. The BMUI and the RMUI cover the polar angles between 34° < 9 < 135° and 
134° < 9 < 171°, respectively. Each chamber has 2 double layers of LST; spatial resolutions of 
1 and 3 mm have been obtained along the direction of the tube axis and perpendicular to it, 
respectively. 

A set of four scintillation counters (C5) in two planes interleaved with a 3 mm Pb foil immedi- 
ately behind the RCAL at approximately Z = —3 m partially surrounds the beampipe. The C5 
counter measures the timing of both beams and also tags events from proton-gas interactions. 
A vetowall (VETO), consisting of two layers of scintillator on either side of an 87 cm thick iron 
wall centred at Z = —7.3 m, was also used to tag and reject off-axis beam particles. 

The ep luminosity was measured from the rate of the Bethe-Heitler process (ep — *> cyp) by 
counting the final state photons in the luminosity monitor. 

1 The ZEUS detector uses a right-handed coordinate system where the Z axis points in the direction of the 
proton beam (forward direction) and the X axis is horizontal, pointing towards the centre of HERA, with the 
nominal Interaction Point (IP) at (0,0,0). 

2 Pseudorapidity is defined as rj — — ln[tan(|)]. 
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4 Kinematics 



The kinematics of the elastic process ep — > e J/ip p are described below. The incoming and 
outgoing electron four-momenta are denoted by k and k', respectively, while the four-momentum 
of the virtual photon is q = k — k'; Q 2 = —q 2 . If P(P') is the incoming (outgoing) proton four- 
momentum, and Pjm is the J/ip four- momentum, then the squared four-momentum transfer t 
is defined as 

t = (P-P') 2 = (q-P^) 2 - 

For Q 2 ~ 0, t ~ — Ptj/^, where prjip is the transverse momentum carried by the J/ip. W is given 
by W 2 = (P + q) 2 . The Lorentz scalar y is defined as P ■ q/P ■ k, and can be approximated as: 

y~{E- pz)j^/2E e , and W 2 = sy ~ AE e E p y ~ 2(£ - pz)j/^E p , 

where E e denotes the electron beam energy, E p is the proton beam energy, y/s is the ep cm. 
energy, and E is the energy and pz is the Z component of the momentum of the J/0- 

5 Trigger and preselection 
5.1 Trigger 

The J/ip was identified from its leptonic decay modes. The momenta of the decay leptons 
from the J/V>'s in the observed kinematic range are low (~ 1.5 GeV). To trigger on these low 
momentum lepton tracks in the high background environment of HERA, where the typical 
background rate from beam-gas interactions exceeds 10 KHz, requires a very selective trigger. 

ZEUS uses a three level trigger scheme^]. The first level trigger (FLT) is built as a deadtime- 
free pipeline. The triggering on the leptons from the decay of the J ftp used the CTD, CAL and 
the muon chambers and is described below. 

Events with low momentum electrons were selected two ways by the FLT: 

• by requiring a minimum energy deposit of 660 MeV in the EMC. To reduce background, 
a total CAL energy deposit greater than 2 GeV or a total energy deposit in the FCAL, 
excluding the region adjacent to the beampipe, greater than 2.5 GeV were required. 
In addition, one to three track segments in the innermost superlayer of the CTD were 
required, 

• by requiring a minimum energy deposit of 464 MeV in REMC and any track segment in 
the innermost CTD superlayer. 

In both of these cases, events were vetoed if the energy deposit in the FCAL region immediately 
surrounding the beam pipe exceeded 3.75 GeV to reduce the background from proton gas 
interactions. 

The analysis of the muon decay mode used the inclusive muon triggers; events with low mo- 
mentum muons were selected by the FLT in one of two ways by requiring: 
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• hits in the RMUI accompanied by an energy deposit of at least 464 MeV in the RCAL, 

• hits in the BMUI accompanied by an energy deposit of at least 464 MeV in the CAL. 

The trigger used only the inner muon chambers to maximize the acceptance. 

All four of these triggers also required at least one CTD track segment pointing towards the IP, 
along with the appropriate VETO and C5 signals to ascertain that the event originated from 
the IP. 

The second level trigger (SLT) reduced the beam related background further by making use of 
the subnanosecond time resolution of the CAL and by requiring that the energy deposits in the 
CAL were in time with the bunch crossing. 

The third level trigger (TLT) ran on a farm of processors; an event was flagged by the TLT as 
a J/ij) candidate if either of the following criteria was satisfied: 

• Electron Mode: a fast electron identification was carried out by using information from the 
CTD and CAL. Electrons were identified by an energy deposit in the EMC of at least 90% 
of the cluster energy, where a cluster is defined as a number of contiguous cells with energy 
deposit. The tracks from the CTD were extrapolated towards the CAL and matched^ with a 
cluster in order to determine the energy deposited in the calorimeter by that track. An event 
was accepted if it had a pair of oppositely charged tracks, each associated with such an electron 
cluster and each with a momentum exceeding 0.5 GeV and px higher than 0.4 GeV. The track 
reconstruction in the TLT used the 'Z-by-timing' information available from the first three axial 
superlayers of the CTD. Additional requirements imposed on the events to remove beam-gas 
interactions were Y^iPzi/^iEi < 0.94 and E^-Ej — pzi) < 100 GeV, where the sums are over all 
calorimeter cells. 

• Muon Mode: for an event with a muon FLT trigger, a match between track segments from 
the CTD and the inner muon chambers was required, with energy deposits compatible with 
those from a minimum ionizing particle in the calorimeter EMC and HAC sections. A track 
with a minimum p? of 1 GeV in the barrel region or a minimum momentum of 1 GeV in the 
rear region was flagged as an inclusive muon trigger and satisfied the TLT requirement. 

5.2 Preselection 

All events satisfying the TLT criteria were reconstructed offline where the more refined in- 
formation from the CTD was used for tracking. The event samples were then divided into 
an electron sample based on kinematics and track-CAL cluster matching, and a muon sample 
based on muon identification by the muon chambers. 

• Electron Mode: events with two oppositely charged tracks, each with px > 0.5 GeV, where 
one of the two had a minimum momentum of 1 GeV, were selected. The tracks, which were 
required to originate from the event vertex, were matched to a CAL cluster as in the TLT, and 
were selected if the invariant e + e~ mass was between 2 and 4 GeV. This data sample contained 
2021 events from a total luminosity of 486 nb _1 . 

3 The distance of closest approach between the extrapolated track and the centre of the cluster was required 
to be less than 30 cm. 
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• Muon Mode: Cosmic ray events, the biggest background in the muon sample, were sub- 
stantially reduced by using the time difference between the upper and the lower halves of 
the calorimeter. Again, events with two oppositely charged tracks, each with px > 0.5 GeV, 
originating from a vertex, were retained. A further reduction in beam gas contamination was 
achieved by requiring Hipzi/^Ei < 0.96, where the sum was over all calorimeter cells. This 
sample contained 456 events from a total luminosity of 490 nb _1 . 



6 Analysis 

6.1 The J /if) signal 

In order to ensure a high quality of track reconstruction, for both data samples only tracks 
within the pseudorapidity range \rj\ < 1.5 were considered. In addition, the total energy in 
the CAL, apart from the energy deposited by the two leptons from the J/^> candidate, was 
required to be less than 1.0 GeV. This criterion was imposed to reject inelastic events. Only 
events where the acceptance was high (the W range between 40 and 140 GeV) were retained. 

• Electron decay mode: a total of 136 events satisfied these criteria with 72 events in the 
invariant mass range 2.85 - 3.25 GeV. The resulting e + e~ invariant mass distribution (M + -) 
is shown in Fig. |2|a. A clear peak is visible at the J/-0 mass. The asymmetric shape is attributed 
to energy loss by the bremsstrahlung process in the material encountered by the decay electrons. 
The same asymmetry was observed in the reconstructed events from the Monte Carlo study 
described below. A maximum likelihood fit to the e + e~ mass spectrum was performed using 
a Gaussian shape convoluted with a bremsstrahlung function to account for the energy loss; 
a second order polynomial in M g + e - was used to describe the background. The processes 
contributing to the background are described in section 6.3. In this fit, the fraction of events 
undergoing bremsstrahlung was constrained to be that determined from the Monte Carlo study 
described below. The mass (3.08 ± 0.01 GeV) and the resolution (38 ± 10 MeV), used as free 
parameters in the fit, were in agreement with the predictions of the Monte Carlo simulation. 
The fit, shown in Fig. |2|a, yields 72 ± 9 events for the J /if) signal. 

• Muon decay mode: the muon data sample still contained some contamination from cosmic 
ray muons. These were identified and removed by applying a collinearity criterion. A total 
of 45 events satisfied the selection criteria in the invariant mass range between 2 and 4 GeV, 
with 35 events in the invariant mass range between 2.85 and 3.25 GeV. Figure ^|b shows the 
invariant mass spectrum: a clear peak at the J/ip mass is observed. The mass spectrum was 
fitted with the combination of a Gaussian and a flat background, shown in Fig. ^|b. The fitted 
mass (3.08 ± 0.01 GeV) and the resolution (61 ± 13 MeV) were again in agreement with the 
Monte Carlo expectations. The fit yielded 32 ± 6 signal events. 



6.2 Monte Carlo simulation and acceptance 

Elastic J/ijj production was simulated with the DIPSI[T5| and the EPJPSIIjTB] generators. The 



model used by DIPSI assumes that the exchanged photon fluctuates into a qq pair which then 
interacts with the proton via the exchange of a pomeron described in terms of a gluon ladder ||. 
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The model is based on a perturbative QCD calculation in the leading log approximation. The 
hard scale in this model is taken to be /i 2 ~ Mjy^/4 ~ 2.5 GeV 2 . The gluon momentum 
density of the proton, used as input in the DIPSI Monte Carlo, was oc x~ 0A . The EPJPSI 
generator assumes pomeron exchange for the elastic J/i/j production. Two different models of 
the pomeron were used. In the first model, the pomeron consists of two gluons, and one of the 
gluons interacts with the cc state into which the photon has fluctuated; the is formed from 
the cc and the remaining gluon of the pomeron. In the second model, the pomeron is without 
a structure, and it interacts with the photon as a whole to form the J/"0. 

Events were generated in the W range between 30 and 200 GeV and between the minimum 
allowed value of Q 2 (~ 10~ 10 GeV 2 ) and 4 GeV 2 . The events were then passed through the 
standard ZEUS detector and trigger simulation programs, and processed with the same recon- 
struction and analysis programs as the data. The distributions of the reconstructed kinematic 
quantities obtained using DIPSI were in good agreement with those from the data. The overall 
acceptance (including the geometric acceptance, detector, trigger and reconstruction efficien- 
cies) was then obtained using DIPSI as the ratio of the number of accepted Monte Carlo events 
to the number generated in the selected kinematic range of W between 40 and 140 GeV. Events 
generated with EPJPSI reproduced some aspects of the data well and were used in the study 
of systematic uncertainties in the acceptance determination. Table 1 shows the acceptances^ 
in various W ranges determined for each decay mode. 

6.3 Background 

Two types of background contributions were considered: the first is a continuum background 
present over the complete mass range, which was already subtracted in obtaining the number 
of signal events from the fit; the second type is J/ifi production through proton dissociation or 
other inelastic production processes. 

• Continuum background: the Bethe-Heitler process produces lepton pairs from photon-photon 
scattering where the electron and the proton each radiates a photon. The invariant mass spec- 
trum of the lepton pair, either e + e~or typically forms a continuum with a maximum at 
low masses. The contribution from this process was obtained by generating Monte Carlo events 
using the generator LPAIR|T7]]. These events were treated the same way as those generated 
by DIPSI. The hatched regions in Fig. |2|a and []b show the LPAIR events, normalized to the 
appropriate ep luminosity, which survived all the selection criteria. In the mass region between 
2.85 and 3.25 GeV, the background from this process in each mode is ~ 8.5%. 

For the electron decay mode, a second source of continuum background was responsible for the 
difference between the quadratic polynomial (dashed curve) and the two photon process shown 
in Fig. |2]a (hatched area). After studies with data and Monte Carlo, this was determined to be 
from pion contamination. One or both of the electron candidates could be misidentified pions; 
the probability of misidentification decreases with increasing momentum. In the mass range 
between 2 and 4 GeV, the background from the pion contamination is comparable to that from 
the two-photon process. 

4 The efficiency of the CAL trigger threshold for the muon decay mode was found to be 90%; this has been 
included in the acceptance values. Also included is an additional factor of 0.9 for the muon chamber efficiency. 
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• the proton dissociation process in reactions like ep — > e J/if> X, where X was undetected, 
can also contribute to the observed J /if) signal, as the outgoing proton is not observed. This 
background was estimated by removing the criterion that the energy deposited in CAL, ex- 
cluding that from the J /if) decay leptons, be less than 1 GeV. A total of four events in the 
electron mode and three events in the muon mode were observed with CAL close to the proton 
direction; these are candidates for J /if) production accompanied by proton dissociation. Monte 
Carlo events simulating this process were generated with PYTHIA[|18|]; the mass distribution 



of the diffractive system X was parametrised as da/dMx 2 ~ Mj"" where n was varied^ be- 
tween 2 and 3. These Monte Carlo events were analysed in the same way as the data. For 
n = 2.5, the proton dissociative contamination in the J/ip signal was estimated to be 17%, by 
comparing the number of events with extra energy in CAL from this Monte Carlo sample with 
the seven events from the data. The fraction of proton dissociative contamination determined 
was (17l| ± 10)%, where the first uncertainty is statistical, obtained from the electron and the 
muon modes together, and the second uncertainty is systematic, observed from the variation 
of n between 2 and 3 in the Monte Carlo. The cross section for the elastic process was thus 
obtained by subtracting 17% from the J /if signal, independent of W. 

• the contribution from the photon-gluon fusion process was determined from Monte Carlo 
studies using HERWIG[plH an d EP JPSlHl6fl . Possible contributions from ^(3685) production 



and processes where the photon undergoes diffractive dissociation e.g., 7 — > 3 /if X, where X 
was not detected, were also considered. The total contribution from these processes was found 
to be negligible and a systematic uncertainty of 3% was assigned to it. 



7 Results 

7.1 ep cross section 

The electropro duct ion cross section, a ep , is calculated as: a ep = ^ x ^: x ^ x ^ where N 
denotes the number of J /if) signal events, A the acceptance, L the integrated luminosity, and Br 
the leptonic branching fraction of J/^p^l, namely (5.99 ± 0.25)% for e + e~and (5.97 ± 0.25)% 
for Subtracting the 17% contribution to the J/if> signal from the proton dissociation 

process, we obtain for the elastic cross section for the process ep — > e J /if p for Q 2 < 4 
GeV 2 and y/s = 296 GeV in the W range between 40 and 140 GeV: a ep = d.6±\i nb, and 
a ep = 6.2tj g nb, from the electron and the muon decay modes, respectively. The errors are 
statistical. 

Figure |3] shows the dcr e p/dpyj^ differential cross section for both decay channels combined after 
the background subtractions. An exponential fit of the form exp(— fcpyj/^) to the distribution 
in the range < Prj/^p — 1 GeV 2 yields a slope b of 3.7 ± 1.0 GeV -2 , while a fit in the range 
< Prj/ip < 0-^5 GeV 2 gives a slope b of 4.5 ± 1.4 GeV -2 , where the statistical and the 
systematic errors have been added in quadrature. The acceptance in these Ptj/^ ranges was 
constant within 15%. As noted earlier, for Q 2 ~ 0, Ptj/^ approximates \t\. Taking into account 
the Q 2 dependence, the slope of the \t\ distribution obtained from the Monte Carlo is ~ 0.5 

5 The CDF Collaboration reported a measurement of n = 2.20 ± 0.03 at y/s = 1800 GeV|l^], in agreement 
with Regge theory predictions |2(J , which is well within the range between 2 and 3 considered. 
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unit higher than the slope obtained from the Ptj/^ distribution, in the Ptj/^ ranges described 
above. 

In order to determine the cross section dependence on W, the data samples were divided into 
two W ranges: 40 < W < 90 GeV and 90 < W < 140 GeV. Table 1 lists the elastic 
cross sections from each W range after subtracting the backgrounds, along with the numbers 
of events and the acceptances. 



Table 1: Acceptance and cross sections 



W range (GeV) 


40 - 90 


90 - 140 


decay mode 


J/ip — > e + e 


J/ip -» fj+n 


J/ip — ► e + e 


J/V' -> 


acceptance 


0.39 


0.11 


0.22 


0.19 


signal events 






16 ± 4 


12±S 


Cop(nb) 


3.9i^±0.7 


a y+i.2+0.9 
' -1.6-1.0 


2.5±0.6±0.4 


9 o +0.9+0.4 
z - z -1.0-0.5 


integrated photon flux 


0.077 


0.033 


er 7p (nb) 


50t? 3 ± 10 


fi1 +16+14 
D1 -20-14 


76±24 ± 13 


fi c + lb+13 
°°-33-14 


(W) (GeV) 


67± 11 


114 ± 9 


a lv (nb) 


52+( 2 ± 10 


71^o ± 12 



7.2 Systematic uncertainties 

The summary of the uncertainties from various sources is reported in Table 2. The polarization 
of the J/0 was not measured and the angular distribution of the decay leptons was varied 
from flat to the form 1 + cos 2 ^*, where 9* is the decay angle of the leptons in the J/0 rest 
frame with respect to the J/-0 lab momentum. This could affect the acceptance by up to 5%. 
No uncertainty was attributed to the W calculation, as the resolution in W is of the order 
of 1%. The variation in acceptance from the modeling of the W dependence was calculated 
using the different Monte Carlo generators described in section 6.2; an uncertainty of 9% was 
assigned to each of the decay modes. The uncertainty in acceptance from the variation in the 
gluon density was included in modeling the W distribution in the Monte Carlo simulations. 
The computation of the muon chamber efficiency added an asymmetric uncertainty as shown 
in Table 2. The uncertainty from the track multiplicity determination was observed to be 5%. 



The energy deposit in CAL was obtained after the uranium noise subtraction |13|. This led to a 
4% uncertainty. The uncertainty due to the calorimeter trigger thresholds for the muon decay 
channel was 10%. For the electron mode, the energy requirement for triggering was far enough 
above the threshold that only a maximum of 5% variation could be observed in the efficiency 
determination. 

An uncertainty of 10% was attributed to the subtraction of the proton dissociation process, and 
a 3% uncertainty was added to account for contributions from other possible inelastic processes 
to the signal, referred to as feed-in from other modes, as discussed in section 6.3. A total 
uncertainty of 17% for the electron mode and ^1% f° r ^ ne muon mode in the electroproduction 
cross section was thus obtained by adding all contributions in quadrature. 
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Table 2: Systematic uncertainties for the ep — ► e J/ift p 
cross section measurement 



uncertainties 


3 /if) — > e + e 


J/V' 


branching fraction 


4 % 


4 % 


luminosity 


3.3 % 


3.3 % 


3 /if) decay angular distribution 


5 % 


5 % 


W dependence 


9 % 


9 % 


muon chamber efficiency 




+5 % -9 % 


track multiplicity 


5 % 


5 % 


<1 GeV energy requirement 


4 % 


4 % 


CAL trigger threshold 


5 % 


10 % 


proton dissociation subtraction 


10 % 


10 % 


feed-in from other modes 


3 % 


3 % 


total 


17 % 


+20 % -21 % 



7.3 Photoproduction cross section 

The photon-proton cross section cr w ^j/^ p is obtained from the corresponding electron-proton 
cross section by using the relation : 



rVmax rQmax 9 / 2\ 

a ep {s) = dy I dQ ■ Q ) • a rp [y, Q 

where 



is the photon flux factor [p3[] , = nil-^, Q ma x = 4 GeV 2 and m e is the electron mass. Since 
the median Q 2 w 10~ 3 GeV 2 is very small, we can neglect the longitudinal contribution and 
the Q 2 dependence of a 7 * p . The cross section is then obtained as the ratio of the measured 
ep cross section and the photon flux factor $ integrated over the Q 2 and y range covered by 
the measurement. This procedure assumes that a lp is independent of y in the range of the 
measurement. As this dependence is not known a priori, the above calculation was repeated 
assuming a rise of a lp (W) proportional to W (= Jsy). An increase of 10% in the resulting cross 
section was found at (W) = 67 GeV and less than 2% at (W) = 114 GeV. These have been 
added in the systematic uncertainty in the photoproduction cross section measurements. The 
cross sections and the integrated photon flux in the different W ranges are summarized in Table 
1. The first uncertainty quoted for the cross sections is statistical and the second is systematic. 
Combining results from the two leptonic decay modes, the measured J/ip photoproduction 
elastic cross sections are : 

Vjp^j/tbp = 52_ + ^ ± 10 nb for (W) = 67 GeV and 
cr 7P ^ J/v , p = 71^ ± 12 nb for (W) = 114 GeV. 

The first error is the weighted error of the uncertainty specific to each decay mode where this 
uncertainty is obtained by combining the statistical error and the systematic error unique to the 
decay mode in quadrature; this error was used to obtain the weighted combined cross section. 
The second error is the systematic error common to both decay modes. 
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8 Discussion 



Figure |] displays a compilation of the J/ip elastic photoproduction cross section measurements [§]. 
The compilation includes only results from the fixed target experiments which measured the 
recoil proton. The elastic J/ip photoproduction measurements by ZEUS are also shown, where 
the elastic cross section has been obtained by subtracting the proton dissociative production 
of the J/ip. A significant rise of the J/ip cross section with the cm. energy is visible. 

Theoretical predictions based on Regge-type and QCD inspired models are also shown in Fig. f|. 
The solid line is the prediction of Donnachie and Landshoff, normalized to lower energy data, 
using a supercritical pomeron||24||. It predicts a slower rise in cross section with W than is 
observed. This is also true of other VDM type models^, ||. The essential point of the QCD 
inspired models is that the cross section is proportional to the square of the gluon density. 
At the HERA energy ranges very low x values (from ~ 5 x 10~ 4 to 5 x 10~ 3 in the present 
analysis) contribute to J/xjj production. The shaded band in Fig. [| is the prediction of the 
Ryskin modelQ, obtained with DIPSI (with a s = 0.25), using the upper and lower limits of 
the leading order (LO) gluon momentum density as extracted by the ZEUS experiment [ 25 ] from 
the scaling violation of F 2 at Q 2 = 7 GeV 2 and evolved back to Q 2 = 2.5 GeV 2 , the scale used 
in the model. The energy behaviour shown by the shaded band is in accord with the data. A 



recent modification by Ryskin [p^j , which uses a calculation based on the present model|| and 
a gluon distribution with some saturation effects is also in agreement with our measurements; 
so is the calculation of Nemchik, Nikolaev and Zakharov|6[], which uses the dipole cross section 
solution of the generalized BFKL ^% equation. 



9 Conclusions 

We have measured the elastic photoproduction cross section of J/ip in ep interactions at 7p 
cm. energies between 40 and 90 GeV ((W) = 67 GeV) and 90 and 140 GeV ((W) = 114 
GeV). The J/^ was detected in its leptonic decay modes (e + e~ or in events where the 

scattered electron and proton were not observed. The elastic photoproduction cross sections, 
obtained from the combined electron and muon decay modes, are : 

o>^j/V>p = 5 2 -il ± 10 nb for (W) = 67 GeV and 
cT^j/^p = 71±H ± 12 nb for (W) = 114 GeV. 

The observed rise in the cross section compared to the lower energy measurements is not ade- 
quately described by Regge-type models and is better represented by perturbative calculations 
if a rise in the gluon density in the proton at low x is assumed. 
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Figure 1: A schematic diagram of elastic 3 ftp production according to (a)VDM with a pomeron 
exchange and (b) QCD-inspired models with the exchange of a gluon ladder. 
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Figure 2: (a) The reconstructed e + e~ invariant mass spectrum for 40 < W < 140 GeV. The 
solid circles represent the data; the solid line indicates a fit to the data with the convolution of 
a Gaussian and a bremsstrahlung function; the dotted line represents a quadratic polynomial 
parametrising the background; the shaded area is the expected contribution from the two- 
photon background, (b) The reconstructed fi + [i~ invariant mass spectrum for 40 < W < 140 
GeV. The solid circles represent the data; the solid line indicates a fit with a Gaussian and a 
flat background. The shaded area shows the flat background expected from the two-photon 
process. 
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Figure 3: the differential cross section da ep /dp 2 T j ,^ for e p — > e J/ip p from both decay 

channels combined for 40 < W < 140 GeV. The continuous line indicates the result of a fit 
with the function exp(— bp^-ii i) with b = 3.7±1.0 GeV -2 . Only statistical errors are shown. 
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Figure 4: A compilation of J ftp elastic cross sections in photoproduction. The solid squares 
represent the measurements from the ZEUS experiment at (W) values of 67 GeV and 114 
GeV. The shaded region represents the prediction of the Ryskin model [5] using the upper and 
lower limits of the gluon momentum density as extracted by the ZEUS experiment at Q 2 = 7 
GeV 2 [25] in LO and evolved back to 2.5 GeV 2 . The solid line is a VDM-like prediction [24]. 
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